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Abstract 

C n 

coefficient of pitching moment 

This paper describes a generic, state-of-the-art, high- 

Cn 

coefficient of yawing moment 

performance aircraft model, including detailed, full- 

Cy 

coefficient of sideforce 

envelope, nonlinear aerodynamics, and full-envelope thrust 
and first-order engine response data. While this model was 
primarily developed for the AIAA Controls Design Chal- 

C 

D 

mean aerodynamic chord, ft 
drag force, lb 

lenge, the availability of such a model provides a com- 

det 

determinant of matrix 

mon focus for research in aeronautical control theory and 

E. 

specific energy, ft 

methodology. This paper also describes an implementation 
of this model using the FORTRAN computer language, as- 

fpa 

flightpath acceleration, g 

sociated routines furnished with the aircraft model, and tech- 

G 

actuator transfer function 

niques for interfacing these routines to external procedures. 

g 

acceleration due to gravity, ft/sec 2 

Figures showing vehicle geometry, surfaces, and sign con- 
ventions are included. 

h 

h 

altitude, ft 

vertical acceleration, ft/sec 2 

Nomenclature 

I 

aircraft inertia tensor, slug-ft 2 

A axial force, lb 

I* 

rotational inertia of the engine, slug-ft 2 

a speed of sound in air, ft/sec 

U 

x-body axis moment of inertia, slug-ft 2 

a B normal acceleration, g 

hy 

x-y body axis product of inertia, slug-ft 2 

Qmj i -body axis accelerometer output, accel- 

h. 

x-2 body axis product of inertia, slug-ft 2 

erometer at center of gravity, g 

Iy 

y-body axis moment of inertia, slug-ft 2 

a*y j/-body axis accelerometer output, accel- 

lyy 

y-2 body axis product of inertia, slug-ft 2 

erometer at center of gravity, g 
a„ z-body axis accelerometer output, accel- 

Iy 

2-body axis moment of inertia, slug-ft 2 

erometer at center of gravity, g 

K 

aerodynamic constant 

a x acceleration along the x-body axis, g 

L 

total body axis aerodynamic rolling moment, ft-lb; 
or total aerodynamic lift, lb 

a v acceleration along the y-body axis, g 

l 

generalized length, ft 

o, acceleration along the 2-body axis, g 

b wingspan, ft 

M 

Mach number; or total body axis aerodynamic 
pitching moment, ft-lb 

C force or moment coefficient 

m 

aircraft total mass, slug 

Cd coefficient of drag 

m. 

mass of engine 

Cl coefficient of lift 

Ci coefficient of rolling moment 

N 

normal force, lb; or total body axis aerodynamic 
yawing moment, ft-lb 

’Electrical engineer, member AIAA. 

n 

load factor 

This paper is a work of the U.S. Government and is not subject to copy- 
right protection in the United States. 

P. 

specific power, ft/sec 



p 

roll rate, rad/sec; or pressure, lb/ft 2 

e 

pitch angle, rad 

Pa 

ambient pressure, lb/ft 2 

p 

coefficient of viscosity 

Pi 

total pressure, lb/ft 2 

p 

density of air, slug/ft 

q 

pitch rate, rad/sec 

ZL 

total body axis rolling moment, ft-lb 

q 

dynamic pressure, lb/ft 2 

ZM 

total body axis pitching moment, ft-lb 

9c 

impact pressure, lb/ft 2 

Z N 

total body axis yawing moment, ft-lb 

Be 

Reynolds number 

4> 

toll angle, rad 

Re' 

Reynolds number per unit length, ft -1 

4>l 

tilt *ngk» of acceleration normal to the flight- 
path from the vertical plane, rad 

r 

yaw rate, rad/sec 

i> 

heading angle, rad 

s 

s 

wing planform area, ft 
complex frequency 

U) 

total rotational velocity of the vehicle 

T 

ambient temperature, K; or total angular 

Superscripts 


momentum slug-ft 2 /sec 2 ; or thrust, lb 

■ 

derivative with respect to time 

T t 

total temperature, K 

Subscripts 

u 

velocity in x-axis direction, ft/sec 

A 

aileron 

V 

total velocity, ft/sec 

or 

aerodynamic reference point 

V c 

calibrated airspeed, kn 

D 

differential stabilator 

v. 

equivalent airspeed, kn 

D 

total drag 

V 

velocity in y-axis direction, ft/sec 

E 

engine 

W 

vehicle weight, lb 

H 

symmetric stabilator 

w 

velocity in z-axis direction, ft/sec 

h 

altitude 

X 

total force along the x-body axis, lb 


measurement not at aerodynamic reference 

X T 

Y 

thrust along the x-body axis, lb 
sideforce, lb 

1 

L 

l 

total lift 
rolling moment 

Y t 

Z 

thrust along the y- body axis, lb 
total force along the z-body axis, lb 

M 

m 

Mach number 
pitching moment 

Z*p 

thrust along the z-body axis, lb 

n 

yawing moment 

a 

angle of attack, rad 

o 

offset from center of gravity 

fi 

angle of sideslip, rad 

V 

roll rate 

Ar 

flightpath angle, rad 

displacement of aerodynamic reference 
point from center of gravity 

r 

9 

R 

pitch rate 
rudder 

Ax 

displacement from center of gravity 

r 

yaw rate 


along x-body axis, ft 

3 

stability axis 

Ay 

displacement from center of gravity 

X 

along the x-body axis 


along y-body axis, ft 

Y 

sideforce 

Az 

displacement from center of gravity 

y 

along the y-body axis 


along z-body axis, ft 

z 

along the x-body axis 

6a 

6d 

6h 

6r 

6jk 

differential aileron command 
differential stabilator command 
symmetric stabilator command 
directional command 
Kronecker delta 

0 

standard day, sea level conditions 
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Introduction 

Background and History 

This paper describes the structure and implementation of 
a high-performance aircraft model. Hie model was devel- 
oped for the AIAA Controls Design Challenge, but is in- 
tended to be useful for a variety of controls and guidance 
applications. Model definition and implementation are cov- 
ered in separate sections. 

The AIAA Controls Design Challenge provides an op- 
portunity for participants to apply control system design 
methodologies to a realistic, nonlinear aircraft model. Any 
rfgrign which performs the control task is acceptable, but in- 
novative or unusual approaches have been encouraged. The 
challenge is a two-year competition. Control designs are 
judged according to their ability to control the model during 
a level acceleration and 3-g turn maneuver at four specified 
flight conditions. 

The model integrates several components. Existing 
pieces were used whenever possible and modified to fa- 
cilitate integration. The result is that most of the model 
implementation is based directly on proven and reliable 
components, although the resulting model is not completely 
representative of any particular aircraft. Because of the 
niiTftri history of the model, some aspects may appear to re- 
semble actual aircraft. The user is warned against making 
any assumptions based on these appearances. 

Model Characteristics 

The model is a collection of modules, each perform- 
ing a specific function. The primary modules are the 
aircraft actuator and surface command inputs, aircraft 
mass and geometry modeling, the equations of motion, 
the atmospheric model, the aerodynamics, the propulsion 


system, and the observation variable modeling. Each ma- 
jor module is described in the following sections. Figure 1 
shows how the modules would be connected together 
with user synthesized control laws to form a complete 
system model. 

Aircraft Description 

The aircraft modeled is a high-performance, supersonic 
vehicle representative of current-day fighters. It is powered 
by two afterburning turbofan engines, each capable of pro- 
ducing approximately 32,000 lb of thrust. A three-view of 
the aircraft is shown in Fig. 2, including control surfaces and 
locations. 



Fig. 2 Three view of aircraft and control surfaces. 

The operational envelope for this vehicle, for trimmed, 
straight-and-levd, 1 -g flight is shown in Fig. 3 for the spec- 
ified weight of 45,000 lb. The envelope includes a maxi- 
mum Mach number (AT) of 2.3 and an altitude limit in the 
50,000- to 60,000-ft range. Mass and geometry parameters 
are given in Table 1 . 

The aircraft primary flight-control surfaces consist of hor- 
izontal stabilators which are capable of symmetric or dif- 
ferential movement, conventional ailerons, and a single 



Fig. 1 Modular structure of the model. 
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Table 1. Mass and geometry characteristics. 


Parameter 

Symbol 

Value 

Wing area 

5 

608.0 ft* 

Wing span 

b 

42.8 ft 

Mean aerodynamic chord 

c 

15.95 ft 

Vehicle weight 
Moments of inertia 

w 

45,000.01b 

(Roll) 

U 

28,700.0 slug/ft 2 

(Pitch) 

I, 

165,100.0 slug/ft 2 

(Yaw) 

i. 

187,900.0 slug/ft 2 

Products 


-520.0 slug/ft 2 


i*. 

0.0 slug/ft 2 


V 

0.0 slug/ft 2 


vertical rudder. The individual surface position limits, rate 
limits, and sign conventions for positive deflection are de- 
tailed in Table 2. The equations in Table 3 define the individ- 
ual surface deflections in terms of command inputs. There 
are a total of five actuators; two aileron, two stabilalor, and 
one rudder. The model includes identical actuators for all 
surfaces. These actuators are rate limited at 24 deg/sec and 
have a first-order response modeled by 


Tfcble 3. Surface deflection definition equations. 


Surface 

Deflection definition 

Rate limit 
deg/sec 


6 A -2 

24 


- 6 a-+2 

24 


( 2Sh — fo) + 2 

24 


(2«h + «d)-h2 

24 

«r 

Sr 

24 


The equations defining the aerodynamic model provide 
nondimensional force and moment coefficients. The lon- 
gitudinal para me ters are in the stability axis system; the 
lateral-directional parameters are given with respect to the 
body axis system. 

The equations used for this model are given in the 
following 

^IrTAM = 4* ^ A* 

= ^Nabc ^ ^ 

Gm*a + Ci MAaK ^No) 

Cd,m = Cd + 


G(s) = 


20 

5+20 


A block dia gram of the actuator model is shown in Fig. 4. 
The command inputs to the aileron and stabilalor surfaces 
are differential and symmetric commands, which are sep- 
arated into inputs to each of the surface actuator models. 
The resulting surface positions are then recombined to ob- 
tain the command response. This is shown for the stabila- 
tors in Fig. 5. Because of the nonlinearities in the stabilator 
command path, the commands will interact in ways which 
are not easily predicted. This interaction is shown in Figs. 6 
and 7. 

Aerodynamic Model 


Cifemv = CyBAMC + Gy, A + Cp+Q fo ^C W|l 

+ Ce tA 6a + Ci h 
-(Ce F P+ Ct 9 r) 


Ct mr = Gg basic + Ge fA 6a + Ct ^ 6d — A Ct i% K^ t 
b 


(2V) 

C% body = Gr basic + G» #A ^A + Gb <d 6d + AG», t K$i u 


( 2V ) 


(C*,p + Cn T r) 


The terms in the equations containing C, A C, A N, or K 
are outputs from the function generation routines, and are 
either calculated by linear interpolation of tabular data or by 
direct calculation. The source of the functional coefficients 
is shown in Table 4. 


The aerodynamics are modeled for the full vehicle en- 
velope using multidimensional tables and linear interpola- 
tion to form nonlinear function generators. In general, these 
aerodynamic quantities are functions of M and some com- 
bination of angle of attack (a), angle of sideslip (/?), and 
symmetric stabilator deflection. 


Propulsion System Model 

The propulsion system model consists of two distinct 
engine models. The engines are similar, but not identi- 
cal; the thrust produced for identical throttle settings is not 


Thble 2. Command input limits and sign conventions. 


Command name 

Symbol 

Limits, deg 

Positive sign convention 

Aileron 

6a 

±20 

Left trailing edge down 

Symmetric stabilator 

6h 

+ 15/— 25 

Trailing edge down 

Differential stabilator 

6d 

±20 

Left trailing edge down 

Rudder 

6r 

±30 

Trailing edge left 
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Mach 810468 

Fig. 3 Vehicle operational envelope at 45,000 lb. 
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Fig. 4 Block diagram of actuator model. 



Fig. 5 Stabilator command path. 
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Fig. 6 Effect of Sh on response. 
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5 d = 0 ° 
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Ifeble4. Source of aerodynamic coefficients. 


Coefficient 

Source 

Independent 

variables 


Table 

M,a,&n 

ACx., 

Table 

M 

CmmAMK 

Table 

M, a, 6h 

AC-.. 

Table 

M 

On, 

Table 

A/, a 

Cm* 

Thble 

M,at 

A No 

Table 

U 

C D (« < 32) 

Table 

Cb*A *C » ^ 

(32 < a < 40) 

Table 

i Of 

(a > 40) 

Calc 

C'l/mAttc i a 

A Co* 

Calc 

h 

ACn ir ... 

Table 

M. throttle 

CtlAK 

Table 


AC t| , 

Table 

M, o,6r 


Table 

M 

Ci, 

Table 

M,a 

Cl, 

Table 

M,a 

Cn ««ac 

Table 

M, a,f) 

c», A 

Table 

M,a 


Table 

M,a 

AC.,, 

Table 

M, a,P 


Table 

M, a 

C. 

Table 

M, a 

Cn, 

Table 

Af, a 

CfiABC 

Table 

M,o,/3 

C» A 

Table 

M,a 

c^ 

Table 

M,ot 

Wo 

AC„, 

Table 

A/,o,6r 

Kh. 

Table 

M 


symmetrical. Each engine thrust vector is aligned with the 
body axis, and acts at a point located 10 ft behind the vehicle 
center of gravity and 4 ft laterally from the centerline. The 
thrust produced by each engine is a function of altitude (k)> 
M, and throttle setting. Each engine is modeled as a non- 
linear system having two separate sections; a core engine 
and an afterburner (augmentor) section with associated se- 
quencing logic. 

Throttle position inputs to the engine model are in de- 
grees, with a mini mum position of 20° and a maximum of 
127°. The core section responds to throttle inputs up to 83°. 
The afterburner section begins to respond at a throttle po- 
sition of 91°. The core model has first-order dynamics and 
rate limiting to model spool-up effects. A block diagram 
of the core model dynamics is shown in Fig. 8. The after- 
burner has a rate limiter and sequencing logic to model the 
fuel pump and pressure regulator effects. A block diagram 
of the afterburner model dynamics is shown in Fig. 9. 


Observation Model 

The observation variables provided by this model repre- 
sent a broad class of parameters useful for vehicle analysis 
and control design problems. These variables include the 
state, time derivatives of state, and control variables. Air- 
data parameters, accelerations, flightpath terms, and other 
mi scel l aneo us parameters are also included. The equations 
used to cal culate those parameters are derived from a num- 
ber of sources (Clancy. 1975; Dommasch et al., 1967; Etldn, 
1972; Gainer and Hoffman, 1972; Gracy, 1980). Implicit 
in many of these observation equations is an atmospheric 
mo del The atmospheric model is derived from the U.S. 
Standard Atmosphere (1962). 

Three body axis angular rates and three translational ac- 
celerations are available as observation variables. These in- 
clude the x-body axis rate (u), the y-body axis rate (v), and 
the 2 -body axis rate (w). The time derivatives of these quan- 
tities, u, ii, an< l w are also included. The equations defining 
these quantities are 

y cos a cos £ 


v = V sin fi 


w = V sin a cos f) 

( Xt — pm sin 0 — Dcos a + L sin a 

• " 

m 

Yt + pm cos flsin 4> + Y 

m 

^ Zt + gm cos 0 cos ^ — D sin a — Lcos 
+qu — pv 



^ + pw — ru 


^ + rti — qw 


The vehicle body axis accelerations constitute the set of 
observation variables that, except for state variables them- 
selves, are most commonly used in aircraft control analysis 
and design problems. These accelerations are measured in 
g units and are derived directly from the body axis forces 
defined in the previous section for translational acceleration. 
The equations used for the body axis acceleration a z , a y , and 
a, are 

a* = (Xt- Dcosa+ Lsin o-ymsin 6)/gom 
a v = (Y T + Y + gmcos 0sin 4>)/gom 

a, = ( Z T - D sia a - Lcos a +gmcos0 cos <(>)/ go m 

where subscript 0 denotes standard day, sea level conditions . 
The equations for the outputs of the body axis accelerom- 
eters (denoted by subscript n) that are at vehicle center 
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Fig. 8 First order engine core dynamics. 
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Fig. 9 Engine afterburner model dynamics. 


of gravity are 

o„ = (Xt - D cos a + L sin a) /go m 
= ( Y T +Y)/gom 

a« = ( Zt - D sin a - Lcos a)/gom 

Ofl = °1U 


relative to the vehicle center of gravity. Included in the set 
of acceleration equations is load factor* n = L/W , where L 
is the total aerodynamic lift and W is the vehicle weight In- 
cluded in the observation variables are the flightpath -related 
parameters, including flightpath angle (7)* flightpath accel- 
eration (fjpa), vertical acceleration (h), flightpath angle rate 
O7)* and scaled altitude rate (h/513). The equations used to 
determine these quantities are 


For orthogonal accelerometers that are aligned with the ve- 
hicle body axes but are not at vehicle center of gravity (de- 
noted by subscript ,t), the following equations apply 

a«x.< = a m - [( q 2 + r 2 )x x - (pq - r) y* 

-(pr+ q)z x ] /go 

a Vti = a^+ [(p«j + f)x„ - (p 2 + r 2 )y„ 

+ (qr-p)z y ] /go 

o«,» = + [(pr — + (qr + p)y, 

~(q 2 + P 2 )z>] /90 


On ,» 


— °nr,» 


7 = sin -1 



fpa = — 

90 

h = a x sin 6 — o„ sin <f> cos 0 — a t cos </> cos 6 
Vh - hV 


7 = 


vVv 2 - /» 2 


Two energy related terms are included in the observation 
variables; specific energy (E,) and specific power (P 4 ), de- 
fined as 

V 2 

E -‘ h * 2? 

9 


where the subscripts x, y, and z refer to the x-, y-, and 
a -body axes, respectively, and the symbols x, y, and z re- 
fer to the x-, y-, and 2 -body axis locations of the sensors 
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The set of observation variables available also includes four 
force parameters; total aerodynamic lift (L), total aerody- 
namic drag (D), total aerodynamic normal force (N), and 
total aerodynamic axial force 04). These quantities are de- 
fined as 

L — qSClm a 
D = qSCDra, 

N = L cos a + D sin a 
A = — L sin a + D cos a 

where Cdwom and C Lmt are coefficients of drag and lift, 
respectively. 

The parameters having the greatest application to 

aircraft dynamics and control problems are the sensed pa- 
rameters and the reference and scaling parameters. The 
sensed parameters are impact pressure (g e ), ambient or free- 
stream pressure (p 0 ), total pressure (p,), ambient or free- 
stream temperature (T), and total temperature (T)). The se- 
lected reference and scaling parameters are Mach number 
(Af ). dynamic pressure ( 5 ), speed of sound (o), Reynolds 
number (Re), Reynolds number per unit length (He'). and 
the Mach meter calibration ratio < 9 e /p«). These quantities 
are defined as 


computed in knots. The calculations assume that internal 
units are in the English system. The equation used for equiv- 
alent airspeed is 

V,= YlMyfi 

which is derived from the definition of equivalent airspeed, 



where po « 0.002378 slug/ft 3 and V, is converted from feet 
per second to knots. Calibrated airspeed is derived from the 
following definition of impact pressure 

_ | [(10+ " l ] (V ' ^ ° o) 

k 1-2 (£) PO [^-6 5 S&/V,y'] -P0 (Ve>0o) 

For the case where V c < oo, the equation for V c is 

V c = 1479.116^ + 1 . 0 ) - 1.0 

Calibrated airspeed is found using an iterative process for 
the case where V c > ao 


a = 


( 


1.4-^rV 

P0T0 ) 


V c = 582.95174 



7.0(Vc/oo) 2 J 






9 = 


pVl 

2 


9c = 


[( 1 . 0 + 0 . 2 M 2 ) 35 - 10 ] p« (M < 1 . 0 ) 


is executed until the change in V e from one iteration to the 
next is less than 0.001 kn. 

The final set of observation variables provided is a miscel- 
laneous collection of other parameters of interest in analysis 
and design problems. The first group consists of measure- 
ments from sensors not located at the vehicle center of grav- 
ity. These represent angle of attack (a ,), angle of sideslip 
(£,), altitude (h/), and altitude rate (/».,) measurements dis- 
placed from the center of gravity by some x-, y-, and z-body 
axis distances. The equations used to compute these quan- 
tities are 



Pt = Pa + 9c 

T t = r(1.0+0.2Af 2 ) 

where l is length, p pressure, T is ambient of free-stream 
temperature, p the density of the air, and p the coefficient 
of viscosity. Free-stream pressure, free-stream temperature, 
and the coefficient of viscosity are derived from the U.S. 
Standard Atmosphere (1962). 

Included in the airdata calculations are two velocities; 
equivalent airspeed (V e ) and calibrated airspeed (Vc), both 


h,i = h+ x sin 0 — y sin <f> cos 0 — z cos 4> cos 6 

h,i = h+ 0(xcos0 + ysin ^sin 0 + zcos ^sin 0) 
—</>(]/ cos $ cos 0 — z sin ^ cos 0) 

The remaining miscellaneous parameters are total angu- 
lar mnmrntnm (T), stability axis roll rate (p,), stability axis 
pitch rate (q,), and stability axis yaw rate (r,), defined as 
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T-J (/,p 2 - 2 IzyPQ - 2 Ixspr + I y q 2 
—2I yt qr + I,r 2 ) 
p s = pcos a + r sin a 

93 = q 

r 9 = — psin a + r cos a 

Equations of Motion and Atmospheric Model 

The nonlinear equations of motion used in this model are 
general six-degree-of-freedom equations representing the 
flight dynamics of a rigid aircraft flying in a stationary at- 
mosphere over a flat, nonrotating Earth. These equations of 
motion were derived by Etkin, and the derivation is detailed 
in Duke, Antoniewicz, and Krambeer. The equations for 
each variable in the state vector are given in the following. 

The following equations for rotational acceleration 
are used 

p= [(XL)Ji + ( X M) h + (ZN)h ~P 2 Ux,h ~ hyh) 

+ pq(I*ih - Ipih - D,h) 

— pr( Iiyh + D v h — lyxh) + ( I y ih — Ixyh) 

— qr(D x I\ - I^h + hih) 

-r 2 (Iy,h -h.h)] / detl 

q = [(X L)h + (ZM)Ia + CZN)h - p 2 Ux,h ~ Ws) 

+ Pq(Ixxh — Iy*I 4 — Dth) 

— pr(.Ixyh + DyU — Iyth) + q 2 (Iyih — Ixyh) 

— qr(D x h — Ixyh + Ixth) 

-r 2 (Iy,I 2 -Ixxh)]/detI 

r= [(XL)J 3 + (XAO/s + (ZN)I 6 -p 2 Uxsh-IxyU) 

+ pq(Ix,h - I y *h - D,h) 

— pr(Ixyh + Dyh — lyxh) + q 2 (Iyth — Ixyh) 

— qr(D x h — I xv Is + Ixth) 

-r 2 (I v ,h-Ix'Is)]/detI 

Where XL, IM and ZN are the aerodynamic total mo- 
ments about the x-, y-. and 2 -body axes, respectively, in- 
cluding power plant induced moments, and 

detl= I x IyI, -llxylxxly, - hl 2 yt ~ ~ Vly 

h = Iyh ~ l] x 

12 = fiyf* + lyxlx* 

J 3 = I zylyx ^ fyfx* 

h = hh - iL 

I5 = + fxyfx* 

h = Ixly - l\y 


D % — It ~ ly 

Dy=I x — I t 
D,= Iy- lx 

The translational acceleration equations used are 

V = [—DcosP + Y sin p + -XV cos a cos p 

+Yt sin f) + Zt sin a cos P 

-my (sin 0 cos a cos 0 — cos 0 sin ^ sin /? 

— cos 0 cos 4> sin a cos fJ)] /m 

a = [— L + Zt cos a — Xr sin a 

+ my(cos 0 cos^cos a + sinflsina)] /Vm cos P 
+5 — tan /J(pcos a + r sin a) 

0 = [Dsin P + Y cos P - Xrcos asm p 

+y r cos P — Zrsin asin P 

+mg( sin 0 cos asin p+ cos 9 sin <j> cos p 

— cos 0 cos ^sin a sin P)] /Vm 
+ psin a — r cos a 

where a, P, 0, and <f> are angles of attack, sideslip, pitch, 
and roll, respectively; Xt , Yt, and Zt are thrust along the 
I-, y-, and z-body axes; and D is drag force, g gravitational 
acc el erati on, L total aerodynamic lift, m total aircraft mass, 

V total velocity, and Y sideforce. 

The equations defining the vehicle attitude rates are 
6 = q cos <f> — r sin <j> 
rp = q sin sec 0 + r cos 4> see 9 
4> = p + q sin <l> tan 9 + r cos <f> tan 6 
where is heading angle. 

The equations defining the Earth-relative velocities are 
h = V(cos0cosasin0-sin0sintf>cos0 

— cos P sin a cos 4> cos 6) 

± = V[ cos /Jcos a cos 6 cos V> 

+ sin 0 (sin 4> sin 0 cos \p — cos 4> sin V>) 

+ cos P sin a(cos 4 > sin 9 cos V> + sin 4> sin VO] 

y = V[ cos p cos a cos 9 sin $ 

+ sin P( sin <j> sin 9 sin V> + cos $ cos 4>) 

+ cos P sin a(cos^sin 9 sin V> - sin 4>cosi>)] 

The atmospheric data model is based on tables from the 
U.S. Standard Atmosphere (1962). This model calculates 
values for speed of sound, acceleration due to gravity, air 
density, viscosity, and ambient static pressure and temper- 
ature. These values are calculated based on altitude. The 


11 



tabular data is organized on evenly spaced breakpoints be- 
tween 0 and 90 km. Linear interpolation is used between 
table values for altitudes in this range; the extreme values 
are used for altitudes outside the range. 

FORTRAN Implementation 

The model is implemented as a program in standard FOR- 
TRAN 77 programming language. The program is intended 
to be portable, and has been successfully compiled on a wide 
range of computer platforms. The routines which imple- 
ment the models are designed in a modular fashion, pat- 
terned after the modules described in the previous section. 
Modules c ommuni cate through named common blocks. 

There are four primary subroutines; each corresponds to a 
pri mar y module in the model. These routines are CCALC, 
ENGINE, INSTRM, and UCNTRL. The CCALC routine 
performs the calculations for the aerodynamic model. The 
UCNTRL routine, supplied by the user, calculates desired 
control surface deflections for the aerodynamic modeling 
routines. The INSTRM routine calculates the instrumen- 
tation models for the observation variables. The ENGINE 
routine contains the propulsion system model. These sub- 
routines and associated common blocks are described in de- 
tail in the following subsections. 

Aircraft Characteristics 

Vehicle mass and geometry information is read from the 
file pun These values are stored in the common 
block DATAIN as follows 

COMMON /DATAIN/ S , B . CBAR, AMSS, 

A IX , AIY , A1Z , AIXZ , 
ADCY, AIYZ, AIXE 

The first three variables in the common block, S, B, and 
CBAR, represent wing area, wingspan, and mean aerody- 
namic chord, respectively. The vehicle mass is represented 
by AMSS. The remaining values are vehicle inertias. 

Equations of Motion and Atmospheric Model Routines 

The routine DERIVC calculates the equations of motion, 
b as**! on the previously listed model equations. The rou- 
tine INTG integrates the equations to calculate the new state 
vector. The integration is performed using a second-order 
Runge-Kutta midpoint algorithm. Atmospheric model cal- 
culations are performed by the ALTFN routine. 

Aerodynamic Model Routines 

The aerodynamic model contains two major subrou- 
tines, ADATTN and CCALC. The subroutine ADATIN is 
u sf/t during ini tialisati on to input the basic formatted aero- 
dynamic data from the file HPAS.dat. The subroutine 
CCALC uses this aerodynamic data, the state variables, and 
the surface positions to determine the aerodynamic coeffi- 
cients. The CCALC subroutine is executed whenever new 
aerodynamic coefficients are required (for example, once 
fjyh frame during a simulation). 


The aerodynamic data are communicated from ADATIN 
to CCALC through named common blocks that occur in 
only these two routines. The interface between CCALC and 
the main program consists of several named common blocks 
that are used to pass state variables, airdata parameters, sur- 
face p o sitions, and force and moment coefficients between 
CCALC and the calling program. 

The main t rans fer of data into the subroutine CCALC 
is through named common blocks. These common blocks 
contain the state variables, airdata parameters, and sur- 
face positions. The transfer of output data from CCALC 
is through a «»"*** common block containing the aerody- 
namic force and moment coefficients. The details of these 
common blocks follow. 

The common block DRVOUT contains the state variables 
and their derivatives with respect to time. The structure of 
this common block is shown in Table 5. The state variables 


Table 5. Names and locations of state 
variables and derivatives 
in DERIVC common block. 


Index 

F 

DF 

Description 

1 

T 

TDOT 

time 

2 

p 

PDOT 

P 

3 

Q 

QDOT 

9 

4 

R 

RDOT 

r 

5 

V 

VDOT 

velocity 

6 

ALP 

ALPDOT 

a 

7 

BTA 

BTADOT 

p 

8 

THA 

THADOT 

e 

9 

PSI 

PSIDOT 


10 

PHI 

PHIDOT 

4> 

11 

H 

HDOT 

h 

12 

X 

XDOT 

x position 

13 

Y 

YDOT 

y position 


are stored first, in the F array, and the derivatives follow 
in the DF array. Each array contains 13 elements. Rou- 
tines which require access to these arrays use an EQUIV- 
ALENCE statement to allow access to the variables by the 
names shown in Table 5, for clarity. For example, the body 
axis rates p, q, and r appear as P, Q, and R, respectively. 
Total velocity is represented by the variable V. and altitude 
by H. 

COMMON /DERIVC/ F(13), DF(13) 

The common Nock SIMOUT contains the main airdata 
parameters required for the function generation subroutine. 
The variables in this common block are 

COMMON /SIMOUT/ AMCH, QBAR. GMA, DEL, 
UB , VB , WB , 

VEAS, VCAS 
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Mach number and dynamic pressure are the first entries in 
the common block, symbolized by AMCH and QBAR, re- 
spectively. The body axis velocities u, v, and w are included 
as UB, VB, and WB, respectively. 

in puts to the actuator models are entered through the 
variables in the CONCOM common block, the structure of 
which is shown in the following 

COMMON /CONCOM/ DAC, DHC, DRC, DDC 

Hie actuator model routine, ACTDYN, calculates the ap- 
propriate inputs to the aerodynamic model, based on the 
command inputs, and stores them in the CONTRL common 
block. The CONTRL common block contains the surface 
position and engine thrust information in the DC array. The 
structure of the common block /CONTRL/ is as follows 

COMMON /CONTRL/ DC (30) 

There are six locations in the DC array used. The defini- 
tion of these is shown in "Ihble 6. Routines which access 
these parameters usually EQUIVALENCE them to the vari- 
able names shown in the table. 


Table 6. Names and locations of command inputs. 


Array 

index 

Command 

variable 

Position 

variable 

Function 

1 

DAC 

DA 


5 

DHC 

DH 


8 

DDC 

DD 

6d 

9 

DRC 

DR 


11 

— 


left engine thrust 

12 

— 

right engine thrust 


The c omm on Mock SIMACC contains the accelerations, 
accelerometer outputs, and normal accelerometer output at 
the center of gravity of the aircraft. 

COMMON /SIMACC/ AX , AY , AZ , 

ANX, ANY, ANZ, 

AN 

The output common block CLCOUT contains the vari- 
ables representing the aerodynamic moment and force 
coefficients 

COMMON /CLCOUT/ CL, CM, CN, CD, CLFT, CY 

The variables CL, CM, and CN are the symbols for the coef- 
ficients of rolling moment (C/), pitching moment (C m ), and 
yawing moment (C«), respectively; these terms are body 
axis coefficients. The stability axis forces are represented 


by CD (coefficient of drag Co). CLFT (coefficient of lift 
C L ), and CY (sideforce coefficient Cy). 

Propulsion System Model Routines 

The top level subroutine ENGINE computes individual 
engine parameters to calculate force, torque, and gyro- 
scopic effects caused by engine offset from the centerline. 
Each time ENGINE is called, it calls the interface routine 
UENGIN, which calls the actual engine model in the rou- 
tine ENGMDL. After ENGMDL computes thrust from each 
engine, ENGINE uses the parameters in the ENGSTF com- 
mon to compute force, torque, and gyroscope effects. The 
control inputs to the engine model are THRSTX(l) for the 
left engine throttle position in degrees, and THRSTX(2) for 
the right engine. The THRSTX array is in the CTPARM 
common block. The output thrust is available in two loca- 
tions; left engine thrust appears as DC(1 1), in the CONTRL 
common block, and as THRUST(l) in the ENGSTF com- 
mon block. Similarly, right engine thrust appears as DC(12) 
and THRUST(2). Various other parameters used for calcu- 
lation of engine effects are passed through common block 
ENGSTF as follows 

COMMON /ENGSTF/ THRUST (4), TLOCAT (4,3), 
XYANGL (4), XZANGL (4), TVANXY (4), TVANXZ (4), 
DXTHRS (4), EIX (4), AMSENG (4), ENGOMG (4) 

The variables in this common Mock correspond to thrust 
(THRUST); the x-, y-, and z-body axis coordinates of the 
point at which thrust acts (TLOCAT); the orientation of 
the thrust vector in the x-y body axis plane (XYANGL) 
in degrees; the orientation of the thrust vector in the x-z 
body axis plane (XZANGL) in degrees; the orientation of 
the thrust vector in the x-y engine axis plane (TVANXY) 
in degrees; the orientation of the thrust vector in the x-z 
engine axis plane (TVANXZ) in degrees; the distance be- 
tween the center-of-gravity of the engine and the thrust point 
(DXTHRS) measured positive in the negative x-engine axis; 
the rotational inertia of the engine (EIX); mass (AMSENG); 
and the rotational velocity of the engine (ENGOMG). 

Although many of the engine arrays are sized to contain 
four elements, only the first two are used. The first element 
corresponds to the left engine, the second to the right engine. 
Engine data is contained in the file ENGINE .dat, stored as 
formatted text. 

Observation Model Routines 

The common block OBSERV allows the user to access 
all the observation variables. This common block contains 
the two vectors OBSVEC(120) and PARAM(120,6). A list 
of the available observations and parameters is given in 
"Bible 7. 
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Table 7. Continued. 


Ifcble 7. Observation variables available 
in the OBSVEC array. 


Location (index) 
in OBSVEC 

Variables 

State variables 

1 

Roll rate 

2 

Pitch rate 

3 

Yaw rate 

4 

Velocity 

5 

Angle of attack 

6 

Angle of sideslip 

7 

Pitch attitude 

8 

Heading 

9 

Roll attitude 

10 

Altitude 

11 

Displacement north 

12 

Displacement east 

Derivatives of state variables 

13 

Roll acceleration 

14 

Pitch acceleration 

15 

Yaw acceleration 

16 

Velocity rate 

17 

Angle-of-attack rate 

18 

Angle-of-sideslip rate 

19 

Pitch attitude rate 

20 

Heading rate 

21 

Roll attitude rate 

22 

Altitude rate 

23 

Velocity north 

24 

Velocity east 

Accelerations 

25 

x-body axis acceleration 

26 

y-body axis acceleration 

27 

x-body axis acceleration 

28 

x-body axis accelerometer at vehicle 


center of gravity 

29 

y-body axis accelerometer at vehicle 


center of gravity 

30 

2 -body axis accelerometer at vehicle 


center of gravity 

31 

Normal acceleration 

32 

x-body axis accelerometer not at 


vehicle center of gravity 

33 

y-body axis accelerometer not 


at vehicle center of gravity 

34 

z-body axis accelerometer not 


at vehicle center of gravity 

99 

Normal accelerometer not 


at vehicle center of gravity 

98 

Load factor 


Location (index) 

in OBSVEC Variables 

Airdata parameters 

91 

Speed of sound 

37 

Reynolds number 

103 

Reynolds number per 


unit length 

35 

Mach number 

36 

Dynamic pressure 

56 

Impact pressure 

55 

Ambient pressure 

57 

Impact/ambient 


pressure ratio 

58 

Total pressure 

59 

Temperature 

60 

Total temperature 

92 

Equivalent airspeed 

93 

Calibrated airspeed 

Flightpath-related parameters 

39 

Flightpath angle 

38 

Flightpath acceleration 

40 

Flightpath angle rate 

43 

Scaled altitude rate 

Energy-related terms 

46 

Specific power 

47 

Specific energy 

Force parameters 

94 

Lift force 

95 

Drag force 

96 

Normal force 

97 

Axial force 

Body axis parameters 

52 

x-body axis velocity 

53 

p-body axis velocity 

54 

z-body axis velocity 

100 

Rate of change of velocity in 


x-body axis 

101 

Rate of change of velocity in 


p-body axis 

102 

Rate of change of velocity in 


z-body axis 

Miscellaneous measurements not 


at vehicle center of gravity 

44 

Angle of attack not at 


vehicle center of gravity 

45 

Angle of sideslip not at 


vehicle center of gravity 

41 

Altitude instrument not at 


vehicle center of gravity 

42 

Altitude rate instrument not 


at vehicle center of gravity 
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Tfcble 7. Concluded. 


Location (index) 
in OBSVEC 

Variables 

Other miscellaneous parameters 

48 

Vehicle total angular 


momentum 

49 

Stability axis roll rate 

50 

Stability axis pitch rate 

51 

Stability axis yaw rate 

Control surface parameters 

61 to 90 

Control surfaces DC(1) 


to DC(30) 

Propulsion parameters 

107 

Throttle position - left engine 

108 

Throttle position - right engine 


The observation vector is produced by the INSTRM rou- 
tine, which contains calculations for instrumentation and ob- 
servation models. The observation vector is in the OBS VEC 
array, in the OBSERV common block as follows 

COMMON /OBSERV/ OBSVEC(120), PARAM(120,6) 

The contents of the observation vector are strictly outputs 
from the simulation; none of the models use these values. 
The stored parameters were selected to provide a variety 
of useful information; however, several unused locations 
(OBSVEC(109) - OBSVEC(120)) are available for the user 
to inrfnrift other observation variables. The user may also 
modify any of the locations without affecting the simulation. 
The PARAM array is used to store parameters which affect 
corresponding measurements in OBSVEC, such as the loca- 
tions of an accelerometer not at the vehicle center of gravity. 

Interfacing and Using the Model as a Simulation 

This section describes the typical method of interfacing to 
and operating the simulation built from the FORTRAN im- 
plementation of the model. The information in this section 
is sufficient to allow a user to interface to the simulation and 
begin using it, without being concerned about the details of 
the model, equations, or the FORTRAN program. 

The routine UCNTRL contains the user control law im- 
plementation. The inputs to this routine may be taken 
from any common block, but typically use the observa- 
tion variables in the OBSVEC array, in the OBSERV com- 
mon block. The OBSVEC array contains many parame- 
ters which would be unmeasurable on an actual aircraft, 
these should not be used as inputs to a practical control sys- 
tem. The outputs are the four control position commands 
in the CONCOM common block, and the two throttle set- 
tings, THRSTX, in the CTPARM common block. A sim- 
ple autopilot is provided in the distribution UCNTRL as an 
example. 


Routines for user supplied input and output are also pro- 
vided. The routine URTIN is called at the start of each frame 
loop, to perform input. The routine URTOUT is called at 
the end of each loop to perform output tasks. These rou- 
tines may perform any function the user wishes. They may 
contain file I/O, interfaces to a network, a windowing sys- 
tem, or shared memory. Using these routines, the user may 
implement displays or recording capability and control the 
simiiiatinn These routines are called in the main loop, so 
their execution will affect execution time significantly. Ex- 
amples of these routines are provided in the distribution. 

A routine is provided for user specified initialization, 
UINIT. This routine is called once, as the last step in 
the initialization sequence. Any initializations may be 
included here. 

The simulation has three mutually exclusive modes of op- 
eration; operate, reset, and hold. The mode selection is per- 
formed by setting one of the variables, OPERAT, RESET, 
or HOLD to true, and resetting the remaining two variables 
to false. These variables are found in the common block 
SSS001. Operate is the normal mode in which the simula- 
tion executes all the input, output, integration, and model 
routines. In hold mode, the simulation only executes the in- 
put and output routines; all states and parameters hold their 
last value. In reset mode, the simulation is reinitialized to 
the user specified initial conditions. The model routines are 
call ed, but no integrations are performed, so simulation time 
does not pass. 

Initial conditions are specified in the array FIC. which 
contains a value for each of the 13 state variables. This array 
is in the common block SETICS. These values are copied 
into the state vector (F) by the ICSET routine whenever the 
simulation is placed into reset mode. Default initial condi- 
tions are provided in the block data initialization, but a new 
set may be used at any time by setting the desired values in 
the FIC array and putting the simulation into reset mode. 

Concluding Remarks 

This paper describes a generic, state-of-the-art, high- 
performance aircraft model, including detailed, full- 
envelope, nonlinear aerodynamics, and full-envelope thrust 
and first-order engine response data. While this model was 
primarily developed for the AIAA Controls Design Chal- 
lenge, the availability of such a model provides a com- 
mon focus for research in aeronautical control theory and 
methodology. This paper also describes an implementation 
of this model using the FORTRAN computer language, as- 
sociated routines furnished with the aircraft model, and tech- 
niques for interfacing these routines to external procedures. 
Figures showing vehicle geometry, surfaces, and sign con- 
ventions are included. 
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